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ABSTRACT 

We present the first Suzaku observation of the new Soft Gamma Repeater SGR 0501+4516, per- 
formed on 2008 August 26, four days after the onset of bursting activity of this new member of the 
magnetar family. The soft X-ray persistent emission was detected with the X-ray Imaging Spec- 
trometer (XIS) at a 0.5-10 keV flux of 3.8 x 10~ n erg s -1 cm -2 , with a spectrum well fitted by 
an absorbed blackbody plus power-law model. The source pulsation was confirmed at a period of 
5.762072 ± 0.000002 s, and 32 X-ray bursts were detected by the XIS, four of which were also de- 
tected at higher energies by the Hard X-ray Detector (HXD). The strongest burst, which occurred 
at 03:16:16.9 (UTC), was so bright that it caused instrumental saturation, but its precursor phase, 
lasting for about 200 ms, was detected successfully over the 0.5-200 keV range, with a fluence of 
~ 2.1 x 10~ 7 ergcm~ 2 and a peak intensity of about 89 Crab. The entire burst fluence is estimated 
to be ~ 50 times higher. The precursor spectrum was very hard, and well modeled by a combination 
of two blackbodies. We discuss the bursting activity and X/7-ray properties of this newly discovered 
Soft Gamma Repeater in comparison with other members of the class. 

Subject headings: pulsar: individual (SGR 0501+4516) — stars: magnetic fields — X-rays: stars 



1. INTRODUCTION 

A new Soft Gamma Repeater (SGR), SGR 0501+4516, 
was discovered on 2008 August 22 by the Swift Burst 
Alert Telescope, thanks to the detection of many short 
bursts (Holland and Sato 2008; Barthelmy et al. 2008). 
Archival ROSAT data showed a faint X-ray source con- 
sistent with the position of this new SGR, probably its 
quiescent counterpart (Kennea & Mangano 2008). 

SGRs are a sub-class of the so called "magnetars", 
believed to be isolated neutron stars with very strong 
magnetic fields (10 14 ~ 15 G: Thompson & Duncan 1995). 
This extreme magnetic field powers their bright emission, 
rather than accretion or rotational power as for most of 
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the X-ray pulsars. Only four confirmed SGRs are known 
to date, all sharing common properties with the other 
magnetars (the Anomalous X-ray Pulsars, aka AXPs) 
such as i) a spin period in a very small range of values 
(2 — 12 s), ii) large period derivatives (10~ 13 — 10~ n s s -1 ), 
iii) bright persistent X-ray emission (10 33-36 ergs" 1 ), 
compared to the rather dim infrared counterparts, iv) 
unpredictable bursting and flaring activity, with a large 
range of energetics (10 37-46 ergs -1 ) and timescales (ms 
to years), and v) transient radio emission (sometimes 
pulsed) connected with their X-ray activity (Woods & 
Thompson 2006; Mereghetti 2008). 

We report in this Letter on a Suzaku Target of Oppor- 
tunity (ToO) observation of SGR 0501+4516, the first 
new SGR discovered in our Galaxy in the last ten years, 
which was performed only four days after the bursting 
activation. 

2. OBSERVATION 

A ToO observation of SGR 0501+4516 was performed 
with Suzaku (Mitsuda et al. 2007), on 2008 August 26 
starting at 00:05 (UT), until August 27 08:25 (UT). The 
X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) 
was operated in the normal mode with 1/4 window op- 
tion, to ensure a time resolution of 2 s. The Hard X-ray 
Detector (HXD; Takahashi et al. 2007) was in the stan- 
dard mode, wherein individual events have a time reso- 
lution of 61 /is, and the four-channel HXD-WAM counts 
are available every 1 s. The target was placed at the 
"XIS nominal" position. 

The XIS and HXD data were both processed with 
the Suzaku pipeline processing ver. 2.2.8.20. Events were 
discarded if they were acquired in the South Atlantic 
Anomaly, or in regions of low cutoff rigidity (< 6 GV for 
XIS and < 8 GV for HXD), or with low Earth elevation 
angles. The net exposures obtained with the XIS and 
the HXD were 43 ks and 55 ks, respectively. 
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Fig. 1. — Left: 2-s binned 0.4—10 keV light-curve of the entire observation, obtained by summing events from all three XIS cameras. The 
largest burst is indicated by a red arrow. Right: A zoom of the ~3 ks of the left panel around the burst. 
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Fig. 2. — XIS results on the persistent emission from 
SGR 0501+4516. Top panel: All XIS spectra of the persistent 
emission (see text for details), simultaneously modeled with an 
absorbed blackbody plus power law. Bottom panel: Pulse profiles 
as a function of the energy band (from top to bottom: 0.5-1, 1-2, 
2-3, 3-5, 5-8, and 8-12 keV). 

3. ANALYSIS AND RESULTS 

3.1. The persistent soft X-ray emission 

The on-source XIS events were extracted from each of 
the 3 XIS cameras, over a region 1'.8 in radius centered 
on the target position. The background events were de- 
rived from a similar region as far away from the source 
as possible. 

We detected SGR 0501+4516 at a 0.4-10 keV count 
rate of 1.74 count s _1 with each XIS front-illuminated 
(FI) CCD (XIS and XIS 3), and 1.82 count s" 1 with 
the back-illuminated (BI) one (XIS 1). The 0.4-10 keV 
XIS light curve presented in Figure Q] reveals 32 short 
(< 2 s) burst episodes, where the count rate per 2 s, 
summed over the 3 cameras, increased by more than 5a 
above the persistent emission (plus background) . Several 
of them were also detected with the HXD ( §3.2|) . 

After eliminating the 32 bursts, we began the timing 
analysis (using Xronos 5.21) by converting all the event 



arrival times to the Solar System barycenter. A power 
spectrum analysis confirmed the reported 5.76 s period- 
icity (Gogiis, et al. 2008) in several harmonics. By fold- 
ing the data at the fundamental period and employing a 
phase-fitting technique (DalPOsso et al. 2003), the best- 
fit period was found to be P = 5.762072 ± 0.000002 s at 
the 14704.0 TJD epoch (all uncertainties are at the 90% 
confidence level). This is consistent with those measured 
with RXTE, Swift, and XMM-Newton (Gogus, et al. 2008; 
Israel et al. 2008a). 

Given the XIS timing resolution of 2 s, we cannot study 
detailed pulse profile sub-structures. However, looking 
at the 0.5-12 keV XIS profiles as a function of energy 
(Figure [2] bottom), we found that below about 5keV one 
sinusoidal function at the fundamental spin frequency fit 
the profile shape well, while at higher energies the second 
harmonic component is needed at a significance of 4cr to 
reproduce the profile. The pulsed fraction (defined as 
the semi-amplitude of the fundamental sinusoidal modu- 
lation divided by the mean background-subtracted count 
rate) was 30(1)% on average, while it varied with energy 
as 21(3), 24(2), 30(2), 43(1), 30(3), and 20(3)%, at 0.4-1, 
1-2, 2-3, 3-5, 5-7, and 7-10 keV, respectively. 

Soft X-ray spectra of the persistent emission were 
studied with XSPEC12. We modeled the background- 
subtracted XIS-FI (XIS plus XIS 3) and XIS-BI spec- 
tra jointly with an absorbed blackbody plus power-law, 
a typical empirical spectral decomposition for magne- 
tars. A multiplicative constant factor was used to ac- 
count for calibration uncertainties between XIS-FI and 
XIS-BI, which were < 3%. As presented in Figure [H we 
found a good fit (xl — 0.96 with 603 degrees of freedom) 
with an hydrogen column density of Nn = (0.89 ±0.08) x 
10 22 cm -2 , a photon-index of T = 2.8 ± 0.1, and a black- 
body temperature and radius of kT = 0.69 ± 0.01 keV 
and 3.2 ± 0.2 dio km, respectively, where d\o is the 
source distance in units of lOkpc. The observed and 
unabsorbed 0.5-10 keV fluxes were 3.8(1) x 10~ n and 
8.7(8) x 10~ n erg s _1 cm~ 2 , respectively. While a good 
fit was also found using a resonant cyclotron scattering 
model (Rea et al. 2008), a two-blackbody model was un- 
successful. The HXD results on the persistent emission 
will be reported elsewhere. 

3.2. The powerful burst and its light curves 

Among the 32 bursts detected by the XIS (§ 13. ip . four 
were also detected by HXD-PIN, and three of them by 
HXD-GSO, all with > 3er significance. As shown in Fig- 
ure [3] with a red arrow, the largest event on 2008 August 
26 at 03:16:16.947 (UTC) was detected by all the Suzaku 
instruments, including the four HXD-WAM. It was also 



SUZAKU OBSERVATION OF SGR 0501+4516 IN OUTBURST 



3 




■100 100 200 300 400 
s -»j_ Time (hie)^ 



03:16 03:17 03:18 03:19 03:20 

Time CUT) 

Fig. 3. — Background-inclusive and dead-time-uncorrected light 
curves of t he la rgest burst, (a) The summed 0.6-10 keV XIS data 
per 2 s (§ ED. (b), (c) HXD-PIN (10-50 keV) and HXD-GSO 
(50-250 keV) data, respectively, with 0.5 s binning, (d) The 28- 
111 keV HXD-WAM data with 1 s binning. Panels (e) and (f) 
are expanded views of (b) and (c), respectively, using 8 ms binning 
(left ordinate), where the 18—1160 keV Konus- Wind light curve (V. 
Palshin, private communication 2008) is superposed in green (right 
ordinate). 



observed by the Konus- Wind instrument (Palmer 2008). 
The HXD-PIN and HXD-GSO light curves are rather 
similar with negligible relative delays (within ±8 ms), 
indicating that the spectrum is not significantly time- 
dependent. 

The HXD-PIN signal intensity at the peak and that 
averaged over the ~ 200 ms duration are ~ 38 and ~ 10 
counts per 8 ms, which translate to a source intensity of 
~ 89 and ~ 23 Crab, respectively. Up to time t ~107ms, 
these light curves are relatively free from instrumental 
dead times (< 6%), which have three components (Taka- 
hashi et al. 2007); processing times in the analog elec- 
tronics (HXD-AE), the limited transfer rate from HXD- 
AE to the digital electronics (HXD-DE), and that from 
HXD-DE to the spacecraft data processor. 

At t ~107ms, the HXD signals suddenly terminated. 
This is an instrumental effect, caused by the second 
and third factors above, including a forced "flush" of 
an HXD-DE output buffer. The number of "lost" HXD 
events can be estimated from various scalar informa- 
tion. Over a 4-s interval including the burst, we found 
that HXD-PIN received ~ 14, 000 photons, of which only 
~ 250, or ~2%, were duly processed onboard and edited 
into the light curve. This fraction is similar in the HXD- 
GSO data, 2.7% (195 processed vs. 7164 received). 

By comparing these Suzaku light curves with the 
Konus- Wind data (green in Fig. [3k and[3f), we found 
that the time history observed by the HXD is a "precur- 
sor" event. When the much stronger main burst began 
at t ~ 138 ms, the HXD data became completely sup- 
pressed, until they returned to normal at t = 740 ms. 
This is consistent with the main burst duration, < 0.5 
s, recorded by Konus- Wind. The entire burst thus con- 
tained ~ 50 times more photons in the hard X-ray band 
than actually detected by the HXD. 

In the XIS data, this burst was split into two readout 



frames of 2 s each. We accumulate the XIS data over 
the 4 s, which contains not only the precursor but also 
the large main peak (missed by the HXD). Around the 
XIS image centroid, the very high burst intensity caused 
strong event pile-up, and telemetry saturation which in 
turn caused some truncation of the XIS frame readout. 
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Fig. 4. — (a) Background-subtracted broad-band spectra of 
the largest burst, fitted by a two blackbody model. Black, red, 
and green specify XIS-FI, HXD-PIN, and HXD-GSO, respectively. 
The ordinate is not divided by exposure, (b) Residuals from the 
fit in panel a. (c)(d) The same as panel b, but when using CompTT 
and CompTT X cyclabs models, respectively, (e) The vF v plot 
corresponding to panel a. The XIS data points are scaled to 6% of 
those in panel (a). 



3.3. The burst spectra 

Figure [3^, shows the spectra of the largest burst, where 
the ordinate is counts per unit energy (not divided by 
exposure). The HXD data were accumulated over a time 
interval of ~ 200 ms around the burst arrival time (see 
Figure |3J). We subtracted the HXD-GSO background 
(22 count s _1 , 50-250 keV) using a 400 s interval be- 
fore and after the burst. The HXD-PIN background 
(0.56 count s -1 , 10-70 keV) was negligible. Thus, the 
emission is detectable up to ~ 200 keV. 

The XIS spectra in Figure [4^ were accumulated from 
XIS-FI over the two frames containing the burst. Un- 
like the HXD data which cover only the precursor, the 
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XIS events must be dominated by the main burst. To 
avoid the pile-up and readout truncation problems, we 
extracted the XIS events using an annulus between 
radii of 2' and 4'. 5, and XIS 3 events between 3' and 4'. 5. 
The two XIS-FI cameras thus yielded 515 burst photons, 
among which pile- up events are estimated as <20%. We 
generated particular XIS ancillary response files that cor- 
rectly reproduce the fraction of photons falling onto these 
limited accumulation regions. 

To evaluate the precursor spectrum, we jointly fitted 
the HXD-PIN (10-70 keV) and the HXD-GSO (50-250 
keV) spectra. The 1-10 keV XIS-FI data were also in- 
corporated, assuming that the precursor has nearly the 
same spectrum as the much brighter main burst. To take 
into account the various uncertainties in the HXD and 
XIS data processing, we left relative model normaliza- 
tions free between them. The interstellar absorption was 
fixed at 0.89 x 10 22 cm -2 as specified by the persistent 
emission. 

As shown in Figure 21 a model comprising two black- 
bodies (typical for magnetars bursts; e.g. Olive et a. 
2004; Feroci et al. 2004; Nakagawa et al. 2007; Israel 
et al. 2008b; Esposito et al. 2008) reproduced the data 
fairly well with Xv = 41.0/35 = 1.17. The deconvolved 
vF v spectra are shown in Figure [4^. The fit yielded 
two temperatures: 3.3lQ4keV and 15.ll 2 'gkeV. Assum- 
ing spherical emission regions, the the cooler and hotter 
blackbodies are estimated to have radii of 8.0^1 <^io km 
and 0.46i°;^ d w km, respectively. The average (over the 
200 ms) and peak 1-200 keV fluxes of the precursor are 
1.0 x 10 -6 and 3.8 x 10 _6 ergs _1 cm~ 2 , respectively, and 
the corresponding fluence is 2.0 x 10~ 7 ergcm~ 2 . The 
total burst fluence is estimated to be ~ 50 times higher. 
The model normalization for the HXD data (PIN and 
GSO together) was - 6% of that of XIS-FI. 

When a thermal Comptonization model, CompTT, is al- 
ternatively used, the fit becomes slightly worse (xv = 
42.8/35= 1.22; Fig.Efc). 

Although the above models are roughly successful for 
the burst spectrum, a hint of negative residuals is present 
at ~ 26keV (Fig. HJd), where the two blackbodies cross 
over. Applying a Gaussian absorption (Gabs) or cy- 
clotron absorption (Cyclabs) factor at ~ 26 keV, the 
fit was improved to \u = 0-89 (Fig. [4ji) . However, its 
addition is significant only at 1.9a level. 

4. DISCUSSION 

SGR 0501+4516 was observed by Suzaku four days 
after the burst activation at a flux level of ~ 3.8 x 



10" 11 ergs" 1 cm 2 (0.5-10 keV). The persistent 0.5-10 
keV spectrum obtained with the XIS is described by a 
kT — 0.69 keV blackbody plus aT~ 2.8 power-law. 

If the identification of SGR 0501+4516 with the 
ROSAT source 2RXP J050107.7+451637 (Kennea & 
Mangano 2008) is correct, the unabsorbed 0.2-10 keV 
flux increased by a factor ~ 15 with respect to the 1992 
value (for a power-law model). Large spectral and flux 
variations have been recorded in other SGRs as well, in 
coincidence with the onset of phases of bursting activity. 
The two most active repeaters, SGR 1806— 20 and, to a 
lesser extent, SGR 1900+14, are known to undergo rather 
gradual variations in luminosity and spectral properties. 
In fact, prior to the Giant Flare of 2004 December 27th 
(Hurley et al. 2005; Palmer et al. 2005), SGR 1806-20 
exhibited a factor ~ 2 flux increase and a decrease in 
T from ~ 2 to ~ 1.6 (Mereghetti et al. 2005). In the 
case of SGR 1900+14 in 2006, the spectral index went 
from ~ 2.4 to ~ 2 and the flux increased by ~ 40% (Is- 
rael et al. 2008b). In contrast, sources which experience 
long stretches of quiescence show much more dramatic 
changes. In May 2008, when SGR 1627-41 reactivated 
after more than 10 yrs, its T switched from ~ 3.3 to 
~ 1.5, and the observed 2-10 keV flux increased by a 
factor of 40 (Esposito et al. 2008). While the luminosity 
increase in SGR 0501+4516 is not dissimilar to that of 
SGR 1627— 41, while in outburst its persistent spectrum 
is much softer. 

Among other short X-ray bursts from SGR 0501+4516, 
wc detected a very powerful one, observed also by Konus- 
Wind (V. Palshin private communication 2008). The 1- 
200 keV precursor spectrum detected with the XIS and 
the HXD was fitted reasonably well by the two blackbody 
model. One possible interpretation of the two tempera- 
tures (~3 keV and ~15 keV) is that they represent the 
photospheres of the ordinary and extraordinary modes 
(Harding & Lai 2006; Israel et al. 2008b), which can 
have different radii and temperatures due to suppression 
of extraordinary-mode scattering and photon-splitting in 
the super-critical magnetic field. 
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